Precise estimation of shell model energy by second order extrapolation method 



Takahiro Mizusaki 
Institute of Natural Sciences, Senshu University, 
Higashimita, Tama, Kawasaki, Kanagawa, 214-8580, Japan 



o 
o 

(N 
O 



> 

in 
o 

<N 
O 
m 
O 



o 

=5 



X 



Masatoshi Imada 

Institute for Solid State Physics, University of Tokyo, Kashiwanoha, Kashiwa, 277-8581, 

(Dated: February 8, 2008) 



Japan 



A second order extrapolation method is presented for shell model calculations, where shell model 
energies of truncated spaces are well described as a function of energy variance by quadratic curves 
and exact shell model energies can be obtained by the extrapolation. This new extrapolation can 
give more precise energy than those of first order extrapolation method. It is also clarified that first 
order extrapolation gives a lower limit of shell model energy. In addition to the energy, we derive 
the second order extrapolation formula for expectation values of other observables. 



PACS numbers: 21.60.Cs 

A pursuit for new numerical algorithms to solve nu- 
clear shell model problems has been one of the most in- 
triguing issues in nuclear structure physics. In the last 
decade, several new methods UHSIllHIHSli] were 
proposed and have been developed. Consequently shell 
model calculation remarkably widens its capability be- 
yond the limitation of conventional approaches such as 
Lanczos shell model diagonalization. Among these new 
methods, several extrapolation techniques U Hi were 
proposed in order to reinforce a restricted power of Lanc- 
zos shell model calculations. We advocated the Lanczos 
shell model diagonalization with the help of the rece ntly 
proposed non-empirical extrapolation technique 0, 0] 
for correlated electron systems and succeeded in enlarg- 
ing the feasibility of large-scale shell model calculations 
by the extrapolation 

In general, interpolation can be reliable if the con- 
sidered quantity behaves smoothly, while extrapolation 
cannot be guaranteed. A reliable extrapolation method 
needs a theory predicting how the considered quantity 
behaves in the extrapolated region. In the recently pro- 
posed extrapolation [i| 0] , it has been shown that, ex- 
act energies defined in the whole space can be extrapo- 
lated from series of energies in truncated Hilbert spaces. 
If these truncated Hilbert spaces give fairly good ap- 
proximations, energies of truncated spaces can be shown 
to behave linearly as a function of the energy variance. 
Therefore exact energies can be estimated by linearly ex- 
trapolated energies of truncated spaces into zero energy 
variance. In condensed matter physics there are good ex- 
amples of this extrapolation in analyses of strongly corre- 
lated electrons on a lattice . In Ref. jl] it has been 
shown that this extrapolation works also when applied to 
the nuclear shell model. However, for detailed spectro- 
scopic studies, precise excitation energies are important, 
and this requires a more precise evaluation of the energy 
eigenvalues. Another problem is, as it will be discussed 
later, that the extrapolated shell model energies are sys- 
tematically lower than the true energy eigenvalues. To 
solve this problem, in the present paper, we successfully 



improve the extrapolation method by including the sec- 
ond order term. Hereafter this extension is called second 
order extrapolation, while the previous one is called first 
order extrapolation. 

In order to present our new method, we consider the 
relation between energy difference and modified energy 
variance as below. We define the difference 5E between 
the lowest energy eigenvalue \H) in a given truncated 
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space and true energy eigenvalue \H) as 
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The modified energy variance AE in the truncated space 
is also defined as 



AE 



H 2 
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Note that this definition of energy variance is different 
from the previous one Q . Throughout the present paper, 
this new definition is used. 

An approximate ground state can be decomposed 
as the sum of the true eigenstate |^o) an( i the rest com- 
ponent \ip r ) as 
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where c, d and c„'s are coefficients, whose normalization 
condition is taken as c 2 + d 2 = 1 and c 2 = 1. The 

IV'o) an d |VVi)' s form an orthonormarized complete set 
spanning the whole space and its energy eigenvalue is 
denoted by E n . By defining the moments Dj (j = 1,2,...) 
as 



c n (En — E ) J 

E 3 
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we can rewrite SE and AE as 



8E = d 2 D 1 E , 



and 



AE = d 2 D 2 E 2 - (d 2 ) 2 (D^o) 2 . 
By eliminating d 2 , we obtain 
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If 4A 2 AE is small enough compared to the unity, we can 
expand SE as 



SE = AA£ + A 3 A£; 2 



(10) 



As the ratio of the moments A embodies the structure 
of the rest component, A can be a function of d 2 which 
is a small parameter. Here we assume that A can be 
expanded by a function of d 2 as 



A = A +A' 1 d 2 + O(<f), 



(11) 



where Aq and A[ are coefficients. This expansion is 
shown to be reasonable in the later discussion. By elim- 
inating d 2 by Eq. (7), A can be rewritten as 



A= Ao + AxAE- 



(12) 



where A\ is a new coefficient for A[ with some additional 
factors. Therefore, we obtain the following second order 
extrapolation formula, 



SE = A AE + (Al + Ai) AE 2 



(13) 



The expansion parameter of Eqs. (12) and (13) is AE/Eq 
but, as Eq is unknown, we present a new extrapolation 
formula by including the Eq depedence in the expansion 
coefficients. We comment that, a dimensionless definition 

of energy variance as AE = AE/ (^H^ leads to the same 
result. 

The first term of Eq. (13) shows the effect of the de- 
crease of the norm of rest component. If we consider 
only the first term, this extrapolation is nothing but the 
previous first order extrapolation 8]. The second term 
shows the second order effects of the decrease of the norm 
of rest component (Aq) and the first order effect of the 
structure change (Ai). The sign of the second order co- 
efficient will be clarified later in the discussion. This new 
method allows us to extrapolate shell model energies in 
the limit AE — > after determining the expansion coef- 
ficients in Eq. (13) by fitting to the shell model energies 
of various truncation spaces. 



In order to investigate the properties of the present 
second order extrapolation method, we utilize numeri- 
cally solvable shell model calculations. One frontier of 
the state-of-the-art large-scale shell model calculations is 
the fp shell. Here we take ground state energies of 48 Cr 
and 49 Cr as fp shell nuclei. For a numerical test, we con- 
sider two cases of different nature. The M-scheme dimen- 
sions are about 2 and 5 Millions for 48,49 Cr, respectively, 
and we can easily diagonalize the Hamiltonian matrices 
by a shell model code, for instance, MSHELL As 
an effective shell model interaction, the KB3 interaction 
is taken. 

We consider two kinds of truncation scheme. The first 
is ©s<t(/7/2) A ~ 40 ~ s ( r ') s where r means the set of the 
/ 5 / 2 , P3/2 and P1/2 orbits and t is the maximum number 
of particles allowed to be excited into the r orbits. This 
truncation scheme is natural due to a considerable shell 
gap between the / 7 / 2 orbit and the others. We call this 
truncation scheme (I). The other is the truncation scheme 

(II): ®s + s'<t ,s'<t(f7/2) A - 40 - s - s '(P3/2) s (r') s ' where r' 
means the set of the / 5 / 2 and P1/2 orbits. Truncation 
spaces are specified by t and to. For 48,49 Cr, the to is set 
to 5, which means up to 5 nucleons in the P3/2 and r' or- 
bits. The parameter t controls the size of the shell model 
space by specifying the maximum number of nucleons in 
the r' orbits. These two sets of truncation scheme grad- 
ually cover the whole shell model space in different ways 
as t increases. 
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FIG. 1: Second order extrapolations of ground state energies 
into zero energy variance; (A) for ground state (0^) of 48 Cr, 
and (C) for ground state ((5/2)]") of 49 Cr. Open circles and 
open squares are energies for the truncated shell model spaces. 
Exact ground state energies are shown by filled diamonds. 
Results of two kinds of truncation scheme are shown for solid 
(I) and dashed (II) curves for (A) and (C). The ratio of the 
monents A is shown as a function of AE for (B) and (D). 
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In Figs. 1 (A) and (C), the energies for the truncated 
shell model spaces with t = 2 ~ 5 for the truncation 
scheme (I) and t = 1 ~ 5 for the truncation scheme (II) 
are shown as a function of the energy variance AE for 
the ground states of 48 Cr and 49 Cr. The solid (I) and 
dashed (II) curves show the second order extrapolations. 
Independently from the nature of the nuclei and from 
the adopted truncation scheme, the shell model energies 
for truncated spaces are well-fitted by quadratic curves 
as indicated by the second order extrapolation formula. 
These extrapolated energies are quite accurate, for exam- 
ples, for 48 Cr, they are -32.91 and -32.97 MeV for (I) and 
(II), respectively, while exact energy is -32.95 MeV. For 
49 Cr, extrapolated energies are -35.59 and -35.60 MeV 
for (I) and (II), respectively, while exact one is -35.59 
MeV. Different truncation schemes can give almost the 
same extrapolation energy. This fact is quite useful for 
confirming the reliability of the extrapolated energy. 
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FIG. 2: First order extrapolations (Li and L2) and second 
order extrapolation (solid curve) are shown for the ground 
state of 49 Cr. 

In Fig. 2, we compare the results of the first order ex- 
trapolation and of the second order extrapolation for the 
ground state of 49 Cr, where the truncation scheme (I) is 
taken, hi and L 2 are the extrapolations from the shell 
model energies for t = 2, 3, 4 and t = 3, 4, 5, respectively. 
L 2 gives higher extrapolation energy than L 1: while it is 
still lower than the exact energy. However, globally the 
shell model energies for truncated spaces can be better 
fitted by a quadratic curve than by straight lines. The 
first order extrapolation is approximately parallel to the 
average tangential line of this quadratic curve around the 
fitted shell model energies. If the quadratic curve is up- 
wardly concave, that is Aq + A\ > 0, the extrapolated 
energy becomes a lower limit of shell model energy and 
approaches the exact energy from the lower side. 

We further go into the details of the second order ex- 
trapolation formula. As we can obtain exact ground state 
wavefunctions for 48,49 Cr, the moments D\ and D2 can 
also be evaluated. In Figs. 1 (B) and (D), the ratio of 
the moments A is plotted as a function of AE. These 
figures show that A can be well described by the linear 
relation of AE. The coefficient Aq is almost the same for 
different truncation schemes, while A\ is different. For 



48 Cr, the moment ratios are A = 0.037 + 0.0018A£; for 
(I) and A = 0.036 + 0.0010A£ for (II). In the second 
order coefficient, A\ is positive for both cases. Moreover, 
A\ is much larger than A\. This means that it is impor- 
tant to take into account the structure change of the rest 
component in second order term. 

As the moments D\ and D2 are positive by definition, 
the ratio of the moments A is positive then Aq is also pos- 
itive. The sign of the second order coeficient of Eq. (13) 
depends on A%. Here we explain the condition A% > 0. 
We define the energy of the rest component relative to 
the true ground state energy, E r = (tjj r \ H \tj) r ) — En. 
The first moment Di is proportional to E r , while the 
second moment D2 is approximately proportional to E%- 
Therefore their ratio A has a 1/E r dependence. If E r 
increases as t increases, A decreases as a function of t, 
which means that A\ > 0. In the following we consider 
two physically relevant but different truncation schemes 
that, however, share a common way for expanding shell 
model spaces as t increases. The shell model orbits are 
divided into higher and lower orbits. As the parameter t 
is the maximum nucleon number in the higher orbits, we 
expand truncated spaces by including basis states with 
larger spherical single particle energies on average. By 
this construnction of the series of truncation spaces, as t 
increases, the rest component is described by basis states 
with larger spherical single particle energies and E r be- 
comes larger. Therefore upwardly concave feature is gen- 
erally expected in energy - energy variance plot. 

Next we consider the expectation value of an observ- 
able O, where a similar second order extrapolation for- 
mula holds. In the same way as the energy, we define the 
difference as 



SO = 10) -10 



that can be rewritten by 

SO ~ d 2 D^E , 
in which D° is defined as 



(14) 



(15) 



jjO ^""^ C n C m (O nrn Ooo^nm) ( l! 



where O nm = (ip. 



En 



ip n ). By eliminating d 2 , SO can 



O 

be expanded by the ratio AE/Eq as 

SO = BAE + (BA 2 ) AE 2 + (17) 

with 

D ? 
D 2 E 

We assume that B can be expanded as 

B = B + BiAE + (19) 
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where Bq and B\ are coefficients. Then, we obtain the 
second order formula for SO as 

60 = B AE + (B Al + Si) AE 2 + .... (20) 

As a further large-scale numerical test, we investigate 
the ground state of 64 Zn, whose M-scheme dimension is 
about 501 Million and its exact shell model calculation 
is quite difficult. However, exact results are reported in 
Ref. [13 . 
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FIG. 3: (A) Second order extrapolations for the energy to 
the zero energy variance for the Of state of 64 Zn. (B, C) 
Extrapolations of the B(E2) (B) and Gamow- Teller (C) sum 
rules for the 0+ state into zero energy variance. Exact re- 
sults are shown by open circles. The effective charges are 
taken as e p = 1.35e and e n — 0.35e. For the oscillator length, 
b = 1.01A 1 / 6 fm is taken. Unquenched Gamow- Teller opera- 
tor are used. The solid and dashed lines show the results of 
truncation schemes (I) and (II), respectively. 

In Fig. 3, the results of extrapolations of energies (A), 
B(E2) (B) and Gamow-Teller (C) sum rules for the Of 



state are shown for 6 Zn. We consider two truncation 
schemes. For the truncation scheme (I), we plot the re- 
sults for t = 1 ~ 4 spaces, whose M-scheme dimension 
is 0.2 ~ 55 Millions. For the truncation scheme (II), 
we plot the results for t = 4 ~ 8 with to — 4, whose 
M-scheme dimension is 0.2 ~ 40 Millions. Shell model 
energies of truncated spaces can be well-fitted by the 
quadratic curves for both the truncation schemes. Ex- 
trapolated energies are -106.38 and -106.34 MeV, for the 
truncation schemes (I) and (II), respectively. The exact 
ground state energy is -106.3 MeV [l5|. The present ex- 
trapolated energies are quite close to the exact one. The 
extrapolated values of the total B(E2) are 1.17 xlO 3 and 
1.18 xlO 3 e 2 /m 4 for the truncation schemes (I) and (II), 
respectively, while the exact value is 1157 e 2 /m 4 [l5j 
with the effective charges as e T = 1.35e and e„ = 0.35e. 
Here we use the second order extrapolation formula for 
the expectation values. The extrapolated values of the 
total B(GT+) with unquenched Gamow-Teller operators 
are 5.3 and 5.5 for the truncation schemes (I) and (II), 
respectively, while its exact value is 5.54. The present 
extrapolations with two truncation schemes can give al- 
most exact values. 

In summary, we have proposed a second order extrap- 
olation method for shell model calculations. Shell model 
energies of the truncated spaces as a function of AE are 
well described by a quadratic curve and the exact shell 
model energy is extrapolated in the limit AE — > 0. To 
derive this second order extrapolation formula, it is quite 
important to consider the second order term of the d 2 
expansion and the structure change of rest components 
in terms of the ratio of the moments A. By numerical 
solvable shell model calculations, we have compared the 
extrapolated energies to the exact one and showed that 
their accuracy is improved by the present second order 
extrapolation and that the first order extrapolation gives 
a lower limit of the shell model energy. In addition, we 
have extended the derivation of second order extrapo- 
lation to any observable. Finally, we have shown that 
different truncation schemes have different behavior in 
extrapolation, but, in principle, extrapolated values must 
be the same. This fact is used to confirm the accuracy 
of extrapolated values. 
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